Microwave-resonance-induced magnetooscillations and vanishing resistance states in 
multisubband two-dimensional electron systems 
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The dc magnetoconductivity of the multisubband two-dimensional electron system formed on the 
liquid helium surface in the presence of resonant microwave irradiation is described, and a new 
mechanism of the negative linear response conductivity is studied using the self-consistent Born 
approximation. Two kinds of scatterers (vapor atoms and capillary wave quanta) are considered. 
Besides a conductivity modulation expected near the points, where the excitation frequency for 
inter-subband transitions is commensurate with the cyclotron frequency, a sign-changing correction 
to the linear conductivity is shown to appear for usual quasi-elastic inter-subband scattering, if 
the collision broadening of Landau levels is much smaller than thermal energy. The decay heating 
of the electron system near the commensurability points leads to magnetooscillations of electron 
temperature, which are shown to increase the importance of the sign-changing correction. The line- 
shape of magnetoconductivity oscillations calculated for wide ranges of temperature and magnetic 
field is in a good accordance with experimental observations. 

PACS numbers: 73.40.-c,73.20.-r,73.25.+i, 78.70.Gq 



I. INTRODUCTION 

The discovery of novel microwave-induced oscillations 
of magnetoresistivity 1 as a function of the magnetic 
field B and so-called zero-resistance states (ZRS)^ has 
sparked a large interest in quantum magnetotransport 
of two-dimensional (2D) electron systems exposed to 
microwave (MW) radiation. The 1/B-periodic oscilla- 
tions were observed for quite arbitrary MW frequencies 
u) larger than the cyclotron frequency uj c . The period 
of these oscillations is governed by the ratio lo/uj c . ZRS 
appear in ultrahigh- mobility GaAs/AlGaAs heterostruc- 
tures as a result of evolution of the minima of the oscil- 
lations with an increase in radiation power. 

Recently, MW-induced magnetooscillations and van- 
ishing of the magnetoconductance a xx were observed 
in the nondegenerate multisubband 2D electron system 
formed on the liquid helium surface 4,5 . These oscillations 
have many striking similarities with those observed in 
semiconductor systems: they are 1/B-periodic, governed 
by the ratio lj/u c , and their minima eventually evolve in 
zero magnetoconductance states nearly at the same val- 
ues oiu)/w c . The important distinction of these new oscil- 
lations is that they are observed only for a MW frequency 
fixed to the resonance condition for excitation of the sec- 
ond surface subband: fno = A2 i (here Ay/ = A; — A;/, 
and A; describes the energy spectrum of surface electron 
states, I = 1,2, ...). 

The ZRS observed in semiconductor systems are 
shownS to be understood as a direct consequence of the 
negative photoconductivity a xx < which can appear 
with an increase in the amplitude of conductivity oscil- 
lations. Regarding the microscopic origin of the oscil- 
lations, the most frequently studied mechanism is based 
on photon-induced impurity scattering within the ground 
subband, when an electron simultaneously is scattered 
off impurities and absorb or emit microwave quanta^. 



This kind of scattering is accompanied by an electron dis- 
placement along the applied dc-electric field whose sign 
depends on the sign of uj c (n — n') + lu (here n = 0, 1, ...). 
Therefore, sometimes this mechanism is termed the "dis- 
placement" mechanism. A different mechanism, called 
the "inelastic" mechanism^, explains conductivity oscil- 
lations as a result of oscillatory changes of the isotropic 
part of the in-plane electron distribution function. 

Both microscopic mechanisms of the negative conduc- 
tivity prosed for semiconductor systems cannot be ap- 
plied for explanation of similar effects observed in the sys- 
tem of surface electrons (SEs) on liquid helium, because 
the MW frequency considered in these theories has no re- 
lation to inter-subband excitation frequencies Ai^/h = 
0Jl,l>- Recently, a new mechanism of negative momentum 
dissipation relevant to experiments with SEs on liquid 
helium was briefly reported 1 ^. It cannot be attributed to 
"displacement" or "inelastic" mechanisms. In this the- 
ory, the origin of magnetooscillations and negative dissi- 
pation is an additional filling of the second surface sub- 
band induced by MW irradiation under the resonance 
condition (tuo = A2,i), which triggers quasi-elastic inter- 
subband electron scattering. The ordinary inter-subband 
scattering, which does not involve photon quanta, is ac- 
companied by electron displacements whose sign depends 
on the sign of uj c (n — n') + U)2,x- Usually, this scattering 
does not lead to any negative contribution to a xx . A 
correction to a xx proportional to lj c (n — n') + U2.1 was 
shown to appear only if N2 > Nie~ A2 > 1 ' Te , where Ni is 
the number of electrons at the corresponding subband, 
and T e is the electron temperature. It is important that 
this correction is also proportional to a large parameter 
equal to the ratio of T e to the collision broadening of 
Landau levels. 

In this work, we perform a systematic theoretical study 
of negative dissipation phenomena in a multisubband 2D 
electron system caused by non-equilibrium filling of ex- 
cited subbands. The magnetotransport theory^ is gener- 
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alized in order to include electron scattering by capillary 
wave quanta (ripplons) which limits SE mobility in ex- 
periments^, where vanishing magnetoconductivity a xx is 
observed. In order to understand the importance of MW 
heating at the vicinity of commensurability points, elec- 
tron energy relaxation is analyzed. A sign-changing cor- 
rection to the energy relaxation rate similar to the sign- 
changing correction to the momentum relaxation rate is 
found for non-equilibrium filling of excited subbands. 



II. GENERAL DEFINITIONS 

Consider a multisubband 2D electron system under 
magnetic field applied perpendicular. The electron en- 
ergy spectrum is described by A; + E n , where e n — 
hoj c (n + I) represents Landau levels (n = 0, 1. 2...). For 
SEs on liquid helium (for review see Ref. Illl ) under a 
weak holding electric field (E± — >• 0), A/ ~ — E R /l 2 , 
where Er is the effective Rydberg energy of SE states, 



E R = 



2m e a 2 B 



, a B 



.A' 



A 



e 2 M 
4(e + l) • 



(1) 



as is the effective Bohr radius, m e is the electron mass, 
and e is the dielectric constant of liquid helium. The 
excitation energy A2.1 is about 6K (liquid 4 He) or 3.2 K 
(liquid 3 He). It increases with the holding electric field 
Ex, which allows also to tune A2,i in resonance with the 
MW frequency. 

Under typical experimental conditions, the electron- 
electron collision rate f e _ e of SEs is much higher than 
the energy and momentum relaxation rates. Therefore, 
the electron distribution as a function of the in-plane 
energy e can be characterized by the effective electron 
temperature, 



f l {e) = N l 2 S-e 



■e/T c 



AZ 



(2) 



where l 2 B = hc/eB, and A is the surface area. Accord- 
ing to the normalization condition J /; (e) Di (e) de = Ni 
[here Di (e) is the density-of-state function for the corre- 
sponding subband], Z\\ = ^„ e _e "/ Te . 

The approach reported here will be formulated in a 
quite general way to be applicable for any weak quasi- 
clastic scattering. As important examples, we shall con- 
sider interactions which are well established for SEs on 
liquid helium. Vapor atoms are described by a free- 
particle energy spectrum = h 2 K 2 /2M with M ^> 
m e . For electron interaction with vapor atoms, it is con- 
ventional to adopt the effective potential approximation 



i41=FW£5>(R e -R ), 



(3) 



where is proportional to the electron-atom scatter- 
ing length 12 . Ripplons represent a sort of 2D phonons, 



and the electron-ripplon interaction Hamiltonian is usu- 
ally written as 



(4) 



where Q q = ■ s /Hq/2puJ q , u> q 



is the ripplon 



spectrum, R e = {z e ,r e }, Sq is the ripplon momentum, 
fe^ q and 6 q are the creation and destruction operators, 
and U q (z e ) is the electron-ripplon coupling^ which has 
a complicated dependence on q. 

For both kinds of SE scattering, the energy exchange 
at a collision is extremely small, which allows to consider 
scattering events as quasi-elastic processes. In the case 
of vapor atoms, it is so because M 3> m e . One-ripplon 
scattering processes are quasi-elastic (huj q <C T) because 
the wave- vector of a ripplon involved is usually restricted 
by the condition ql B < 1. 

For quasi-elastic processes in a 2D electron system un- 
der magnetic field, probabilities of electron scattering are 
usually found in the self-consistent Born approximation 
(SCBA) 13 . Following Ref. [H, we shall express the scat- 
tering probabilities in terms of the level densities at the 
initial and the final states. Then, Landau level densi- 
ties will be broadened according to the SCBA 13 or to the 
cumulant expansion method^, 



Di (e) = - 



A 



2ir 2 l 2 B h 



^2 Im G Un (e) , 



(5) 



where Gi :Tl (e) is the single-electron Green's function. 
The later method is a bit more convenient for analyti- 
cal evaluations because it results in a Gaussian shape of 
level densities 



ImG;^ (e) = 



2ttH 



r. 



1,11 



■ exp 



2(e- 



r 2 



(6) 



Here T^ n coincides with the broadening of Landau levels 
given in the SCBA. For different scattering regimes of 
SEs, equations for T^ n are given in Ref. [llj. We shall 
also take into account an additional increase in T2 : n due 
to inter-subband scattering. 

Effects considered in this work are important only un- 
der the condition r; „ -c T which is fulfilled for SEs on 
liquid helium. Therefore, we shall disregard small correc- 
tions to Z\\ caused by collision broadening because they 
are proportional to Tf n /8T 2 . In other equations, some- 
times we shall keep terms proportional to Fi iH /T e , if they 
provide important physical properties. 

Average scattering probabilities of SEs on liquid he- 
lium and even the effective collision frequency i/ e g can 
be expressed in terms of the dynamical structure factor 
(DSF) of the 2D electron liquid^ S(q,u). This proce- 
dure somehow reminds the theory of thermal neutron (or 
X-ray) scattering by solids, where the scattering cross- 
section is expressed as an integral form of a DSF. With- 
out MW irradiation, most of unusual properties of the 
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quantum magnetotransport of SEs on liquid helium are 
well described by the equilibrium DSF of the 2D elec- 
tron liqui d 11 ^ 16 . A multisubband electron system is ac- 
tually a set of 2D electron systems. Therefore, the single 
factor S (q, ui) is not appropriate for description of inter- 
subband electron scattering. Luckily, for non-interacting 
electrons, we can easily find an extension of S (q, oj) which 
could be used in expressions for average scattering prob- 
abilities of a multisubband system: 

Sl,v M = ^zY<j dee-^Jl n ,(x q ) x 

x ImG^n (e)IniGV,n' (e + fiw), (7) 



where 



x(n,n')}\ 



[max(n, n')] 



T \ n ~ n '\ t~\ 



x q = q-l 2 B /2, and L™{x) are the associated Laguerre 
polynomials. The factor Si^ (q, oS) contains the level den- 
sities at the initial and the final states, and it includes 
averaging over initial in-plane states. At / = V , this fac- 
tor coincides with the DSF of a nondegenerate 2D system 
of non-interacting electrons. Generally, Siy (q,w) is not 
the dynamical structure factor of the whole system, nev- 
ertheless this function is very useful for description of 
dissipative processes in presence of MW irradiation. 

As a useful example, consider the average inter- 
subband scattering rate i>\^v caused by quasi-elastic 
scattering, which is important for obtaining subband oc- 
cupancies n\ — Ni/N e under the MW resonance 17 . Using 
the damping theoretical formulation^ and the SCBA 13 , 
can be represented in the following form 



q 



Xl,V (q)Si,v (q,uu>), 



(8) 



where XU'( = Xl',l) describes electron coupling with scat- 
terers. For SEs on liquid helium, we have two kinds 
of scatterers: ripplons and vapor atoms. Therefore, 
Xi,v = Xu' + Xii" Electron-ripplon scattering gives 



where N„ 



—0 2 2N 



m e T 
ah 3 q 2 



( e WT_!) 1 > 1; and (jj q ) iv 
(l\ U q (z e ) \V). For electron scattering at vapor atoms, 



(9) 



Xi,i> (<?) = Pi,v, v o 
where 



,( Q )„ 



(10) 



B 



Pi, i' 



i,i 



Biu 



■K, 



'I'd 



L z is the height above the liquid surface, ni 31 ^ is the den- 
sity of vapor atoms, and K z is the projection of the vapor 



atom wave-vector. The represents the SE collision 
frequency at vapor atoms for B = 0. 

The generalized factor of a multi-subband 2D electron 
system Siy (q,ui) will be used throughout this work be- 
cause its basic property 



S LU (q,-Lj) = e- hu / T °S v<l (q,u>) 



(11) 



allows us straightforwardly to obtain terms responsible 
for negative dissipation. This property follows from the 
detailed balancing for quasi-elastic processes, w>-yi = 



-A,,//T e 



Vi-^i' , and also directly from the definition of 



Eq. (|7J| . Using Gaussian level shapes of Eq. 
find 



one can 



2ir 1 / 2 h 



Z 



where 21^.,,^ = I*, 



— e en/1 °Il,n;l>,n>(u), 

, J- l,n\l' ,n' 

(12) 



exp 



rf'.n'j and 
Il,n;V,n+m ( w ) = 
fojj — mhio c — Tf n /4T e 



I' ,n+m 



r 2 

1 l,n 

m 



.(13) 



The Si } i> (q,u), as a function of frequency, has sharp 
maxima when u) equals the in-plane excitation energy 
(n! — n) haj c . The parameter r 2 n . ; , , describes broad- 
ening of these maxima. Eqs. (fT2|) and (fT3|) satisfy the 
condition of Eq. (fTTj) . Terms of the order of (IY„/T e ) 
entering the argument of Eq. (fT3| could be omitted, as it 
was done for Z||, because even the linear in Ti jn /T e term 
provides us the necessary condition of Eq. (fTTj) . Anyway, 
our final results will be represented in forms which al- 
low to disregard even the linear in Ti jn /T e term entering 

Il,n;l',n> ( w )- 

Consider the decay rate of the first excited subband 
£*2->-i- Under typical experimental condition, Ti, n; i^ n i is 
much smaller than hio c . Therefore, most of terms en- 
tering Siy (q, cj2,i) are exponentially small and can be 
disregarded. The exceptional terms satisfy the condition 
n' — n = to* (B), where to* (B) = round (w2,i/^c) is an 
integer nearest to u>2,i/w c . In this notation, 



*2-H = ]T 
n=0 



2,n;l,n+m* 



Z\\ K 1 / 2 Y 2 ,n;l,n+r. 



x exp ■ 



ft 2 (^2,1 — m*uj c Y 
V 2 



(14) 



where 



Pl,n;V ,n+m — I Xl,l' Jn,n+m { x q)dx q . 
JO 

For electron scattering at vapor atoms, fti, n ;i',n 



/9j = VqPiji which coincides with Xu> ■ I n the case of 
electron-ripplon scattering, P2,n;i,n+m has a more com- 
plicated expression due to a particular form of U q {z e ) 
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entering the definition of Xii'- The i>2^x(B) is a 1/B- 
periodic function. It has sharp maxima when W2,i/w c 
equals an integer. In the argument of the exponential 
function of Eq. (1141) . we have disregarded terms which 
are small for Ti^ n /T e <C 1. 

Transition rates vi^i' determine subband occupancies 
ni under the MW resonance. At low electron tempera- 
tures, the two-subband model is applicable, and the rate 
equation gives 



Tl2 

ni 



r + D 2 -+i 



(15) 



where r is the stimulated absorption (emission) rate due 
to the MW field, and m + n 2 = 1. Thus, under the MW 
resonance, magnetooscillations of p2->i lead to magne- 
tooscillations of subband occupancies fi\ and n 2 . For 
further analysis, it is important that MW excitation pro- 
vides the condition n 2 > nie _A2 ' 1 / Te , which is the main 
cause of negative momentum dissipation. 



III. MAGNETOCONDUCTIVITY UNDER 
RESONANCE MW IRRADIATION 

Consider now an infinite isotropic multisubband 2D 
electron system under an in-plane dc-electric field, as- 
suming arbitrary occupancies of surface subbands ni in- 
duced by the MW resonance. In the linear transport 
regime, the average friction force acting on electrons due 
to interaction with scatterers F sca t is proportional to the 
average electron velocity V av = (v e ). This relationship 
can be conveniently written as F sca t = — N e m e v c ftV av , 
where the proportionality factor v c s represents an effec- 
tive collision frequency which depends on B and, gener- 
ally, on electron density. The F scat is balanced by the av- 
erage Lorentz force (Ffidd), which yields the usual Drude 
form for the electron conductivity tensor cr^fc, where the 
quasi-classical collision frequency is substituted for 

,/ -11.16 

The effective collision frequency v c f{ can be obtained 
directly from the expression for the average momentum 
gained by scatterers per unite time. Usually, to describe 
momentum relaxation, one have to obtain deviations of 

I 



the in-plane electron distribution function from the sim- 
ple form of Eq. ([2} induced by the dc-electric field. For 
the highly correlated 2D system of SEs on liquid helium 
under magnetic field, this problem was solved in a general 
way, assuming that in the center-of-mass reference frame 
the electron DSF has its equilibrium form (q,uj). In 
the laboratory frame, its frequency argument acquires 
the Doppler shift S (q, w) = (q, w — q ■ V av ) due to 
Galilean invarianc o 11 ' 16 . This approach is similar to the 
description of electron transport by a velocity shifted 
Fermi-function of the kinetic equation method, where 
V av is found from the momentum balance equation. The 
same properties can be ascribed to the generalized fac- 
tori2S M , (q,w). 

Here, we consider a different way, taking into account 
that F scat , as well as the momentum gained by scatter- 
ers, can be evaluated in any inertial reference frame. We 
choose the reference frame fixed to the electron liquid 
center-of-mass, because in it the in-plane distribution 
function of highly correlated electrons (v ee ^> z/ e g) has 
its simplest form of Eq. @, and the generalized factor 
Si ji (q, lj) has it equilibrium form of Eq. ([7]). Then, F scat 
can be considered as the drag due to moving scatterers. 
At the same time, distribution functions of scatterers 
which arc not affected by external fields can be easily 
found according to well-known rules. 

In the electron liquid center-of-mass reference frame, 
the in-plane spectrum of electrons is close to the Lan- 
dau spectrum, because the driving electric field E' = 
E — (l/c)B x V av is nearly zero, at least for uj c v e g. 
At the same time, in this frame the ripplon excitation 
energy is changed to Eq = Hu)q — /iqV av , because the 
gas of ripplons moves as a whole with the drift veloc- 
ity equal to — V av . The same Doppler shift correction 
— fiqV av appears for the energy exchange in the case of 
electron scattering at vapor atoms, even for the limiting 
case M — ¥ oo (impurities which are motionless in the 
laboratory reference frame). In the electron center-of- 
mass reference frame, vapor atoms move with the veloc- 
ity — V av and hit electrons which results in the energy 
exchange — hqV av . 

Describing electron-ripplon scattering probabilities in 
terms of the equilibrium factor Siy {q,uj), as discussed 
above, contributions to the frictional force from creation 
and destruction processes can be found as 



(16) 



r 



It is clear that disregarding the Doppler-shift correction 
— fiqV av in this equation yields zero result. This cor- 
rection enters the ripplon distribution function iV q and 
the frequency argument of the factor Sy. In the linear 
transport regime, the Doppler-shift correction entering 



the ripplon distribution function is unimportant. This 
can be seen directly from Eq. (|16|): setting Eq J — > 
in the frequency argument of gives zero result for 
Fscat- Therefore, in this equation one can substitute 7V q 
for N q defined in Eq. ©. We can also disregard huj q in 
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the frequency argument of Sij'. Then, interchanging the 
running indices of the second term in the square brack- 
ets, and using the basic property of Sy/ (q,u>) given in 
Eq. (fTTj) , Eq. (|T5|) can be represented as 



2m P A 



^2 h< iXi,i' (q) Sim (q, wi,i> + q • V a 



1,1' q 



x (ni - nre- A '.''/ r -e- R * v «/ r ') , (17) 

where Xl,l' (?) = (?) was defined in Eq. fljJJ). This 
equation has the most convenient form for expansion in 

fiq- v av . 

A similar equation for F scat can be found consider- 
ing electron scattering at vapor atoms. Evaluating mo- 
mentum relaxation rate, one can disregard Hwk. k' = 

£k'_k — £ K' w hich represents the energy exchange at a 
collision in the laboratory reference frame. In the center- 
of-mass reference frame, Doppler-shift corrections enter 
the vapor atom distribution function N^) and the fre- 
quency argument of the factor Sim due to the new energy 
exchange at a collision — hqV av . The correction entering 
Nj£) is unimportant because of the normalization condi- 
tion: J2k> ^k' = n a L Z A. Therefore, we have 



N„hv, 



(a) 



.A 



V x 



l,V 



av/ ' 



(18) 



In order to obtain the form of Eq. (|17j) , we represent the 
right side of Eq. (|T5|) as a sum of two identical halves and 
change the running indices in the second half: q — > q 
and £ Then, the basic property of Sim (q,u>) yields 

Eq. (fT?) with xz,i' = xj"/, where x[°! is from Eq. (|TU1) . 

Thus, Eq. (fl7|) is applicable for both scattering mecha- 
nisms. In the general case, xi,l' — Xu> + Xii" The effec- 
tive collision frequency under magnetic field v b q can be 
found expanding Eq. (fP7| in q • V av up to linear terms. 
We shall represent u c s as a sum of two different con- 
tributions: t'cff = i^n + ^A- The normal contribution 
i/N originates from the expansion of the exponential fac- 
tor exp (— h(\ ■ V av /T e ). In turn, ifa can be represented 
as a sum of contributions from intra-subband and inter- 
subband scattering z/n = vn, intra + ^n. inter- The sums 
°f ^N.inter take account of all 1,1'. It is useful to rear- 
range terms with I < I' (Aim < 0) by interchanging the 
running indices I <=^ I', and using the basic property of 



Sim (q,u). Then, we have 
hw 2 



.intra 



.inter 



foJ f 

TY^l™ 1 / x qXl 1 i(Q)Si i i(q,0)dx q , (19) 



47rT f 



£( 

i>i' 



-A.J//T, 



x< / x q xiM (q) Sij' (q,u>ij')dx q . (20) 

The z/n (B) is always positive. In the limiting case of 
a one-subband 2D electron system (hi — Si t \), Eq. (fT9|) 
reproduces the known relationship between the effective 
collision frequency and the electron DSF— . In the paren- 
theses of Eq. (l20l) , the first term is due to scattering from 
I to V, while the second term describes the contribution 
of scattering back from I' to I. It should be noted that the 
forms of Eqs. (fl9)) and (f20f allow to simplify Sim (Qi^im) 
of Eq. (|12p by disregarding small corrections proportional 

to T Ln /T e and (Ti. n /T e ) 2 entering Ii,n;l',n' (w) defined by 
Eq. (fTJ). 

The anomalous contribution to the effective colli- 
sion frequency va (B) can be found from Eq. (1171) ex- 
panding Si t i> (q, wiji + q • V av ) in q-V av , and setting 
exp (— ftq • V av /T e ) — > 1 in the parentheses. In this case, 
to rearrange terms with I < V (Aim < 0)' we shall use 
the property 



(?, -w) = -e r:^ l/ (g jW ) + _ e ^^(g^) 



(21) 



Here 5;;/ (<7,w) = dSi,v (<7i w ) and the last trans- 
formation assumes that r^ n -C T e . Interchanging the 
running indices I I' of terms with Ay < and using 
Eq. (|2Tj) . (5) can be found as 



i>v 



m - ni>e 



-A,,//T e 



< / £<?XM' ( < ?) s 'u' (q^Lv)dx q . 
la 



(22) 



As compared to v^,inter of the normal contribution, here 
the second term in parentheses has the opposite sign. 
Therefore, for usual Boltzmann distribution of subband 
occupancies, isa (B) — 0. The anomalous contribution 
appears only when hi ^ hi'e~^ 1 - 1 ' ^ Te , which occurs under 
the MW resonance condition uj = 0Jn'. 

In the form of Eq. (|2"2")l . it is possible to use a simplified 
expression 



2 

n,n' 
7i,n' 



»/T e 



/' ,n' 



exp 



[Aim - ( n ' - n ) ^c] 2 { 2fr 2 [ujim - (n' - n) u c 



p2 



r 2 



(23) 
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which disregards terms proportional to Yi. n /T e and 
{T Ln /T e f. From Eqs. ([22]) and ((23j) one can see that 
at hi > hi'e~ A '- 1 ' / T " , the sign of va (B) is opposite to the 
sign of loiji — [n 1 — n) lo c . Therefore, va (B) < when the 
magnetic field B is slightly lower the commensurability 
condition A2 i i/Sw c = m (here m is an integer), which 
agrees with the experimental observation for minima of 



For further analysis, it is convenient to introduce 



A 



l.rid' 



XqXl,l> (<?) J'n,n' [x q ) dx q 



(24) 



In the ultra-quantum limit hui c T e , terms with n > 
entering Eq. (|27j) can be omitted, which allows to de- 
scribe magneto-oscillations of va (B) in an analytical 
form. In contrast with oscillations of the normal contri- 
bution in the vicinity of the commensurability con- 
dition, va is an odd function of loi.i/^c — W*- 

Thus, the effective collision frequency f c ff = + 
and magnetoconductivity a xx of SEs are found for any 
given electron temperature. In order to obtain T e as a 
function of the magnetic field, it is necessary to describe 
energy relaxation of SEs for arbitrary subband occupan- 
cies. 



When referring to a particular scattering mechanism, we 

shall use a superscript, \i >n -i',n' = ^fXi',n' + ^%l-,l',n'- 
Consider a two-subband model which is valid at low 
enough electron temperatures. Using the new definitions 
given above, the normal contribution to the effective col- 
lision frequency can be represented as 



n=0 



-e n /T e 



(hu c y 



2V?T e Z|| 



»1- 



A, 



n 2 - 



(25) 



^N, inter 



n 2 + me 



\2,n;l,TC+m* 
■ 2, n; 1,71+771* 



i/T e 



(fee 



-e»/T e 



■ exp 



h 2 (o;2.i — m*uj c ) 

2, n;l, 71+771* 



,(26) 



where m* = round (w 2 ,i/ w c) is the function of B de- 
fined in the previous Section. The j/n, intra and ^n, inter 
have magnetooscillations of two kinds. Oscillations of 
Vn, inter are quite obvious, because quasi-elastic inter- 
subband scattering increases sharply at the commensu- 
rability condition: u>2,i — mu) c . The shape of these peaks 
is symmetrical with respect to the point A.2,i/fui) c = m. 
It is formed by the interplay of the exponential factor, 
having T2 t n;i,n+m* for the broadening parameter, and 
the line-shapes of the subband occupancies. It should 
be noted that at low electron temperatures, ^n, inter is 
exponentially small. The intra-subband scattering con- 
tribution v^ t intra oscillates with 1/B in an indirect way 
because of oscillations in level occupancies hi and h\ in- 
duced by oscillations in the decay rate ^2->i, according 
to Eqs. (fT4|) and (fT5| . These oscillations have also a 
symmetrical shape whose broadening is affected by the 
relation between r and i?2— s-i • 

Magnetooscillations of va {B) have a completely dif- 
ferent shape: 



va 



ri2 — n\e 2,1/ ° 



{hw c ) 



n=0 



p -e n /T e \ 

e / »2,n;l,n+m* 
1 2,n;l,n+m* 



x exp 



fo 2 (^2,1 — m*ui c Y 

2,n;l,n+m* 



2h (w2,i — m*oj c ) 



r 



2,7i; l,n+m* 



(27) 



IV. ENERGY DISSIPATION 

It is instructive to analyze another important example 
of negative dissipation which can be induced by the MW 
resonance. Consider the energy loss rate of a multisub- 
band 2D electron system due to interaction with scat- 
terers. In this case, there are no complications with the 
dc-driving electric field or with the Doppler shifts which 
can be set to zero. This analysis will be important also 
for description of electron heating due to decay of the SE 
state excited by the MW. 

The energy loss rate per an electron due to one-ripplon 
creation and destruction processes can be represented in 
terms of Sj,// (q, oj) quite straightforwardly: 



(u q ) 



l,V 



niSi,i> (q, ui,v - uJ q ) - nie hUq/T 'S u * (q, uj im + oj q ) 

(28) 

Interchanging the running indices in the second 

term, and using the basic property of Si.i' (q,uj) given 
in Eq. (|lip. the terms entering the square brackets can 
be rearranged as 



ni _ ni , e -A hll /T ee -^ q (l/T-l/T e ) 



Sim (q,ui,i> - uj q ) 

Since the processes considered here are quasi-elastic, we 
can expand this equation in frw q and represent W as a 
sum of two different contributions: W = Wn + Wa- The 
normal energy loss rate Wn is proportional to T e — T, 
which is a measure of deviation from the equilibrium, 



W N 



Here 



(30) 



q 1,1' 



~(r) 



2phT e 



(Uq)u> 



This contribution originates from expansion of the expo- 
nential function in huj q (1/T — l/T e ). 
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It is conventional to represent the energy loss as Wn = 

(r) (r) 

— (T e — T) i? N , where is the energy relaxation rate 
of an electron. Rearranging terms with I < I' (A/,p < 0), 
as described in the previous Section, one can find 



~(r) 



e q l 



m e A 



q 1>V 



(31) 

~(r) 

The normal contribution is always positive, which 
means positive dissipation (Wn < 0) regardless of sub- 
band occupancies fii . 

An anomalous contribution Wa appears when ex- 
panding Su' {q,ujiji — <jj q ) of Eq. (|29|) in uj q and setting 
e -hui q (i/T~i/T e ) ^ Q um ty. The rearrangement of terms 
with I < V (A/,// < 0) based on the property of Eq. (j2"Tj) 
yields 



2TT„ 



1>V 



(32) 



Here n; 



-A.,,/r e 



represents an additional mea- 



sure of deviation from the equilibrium induced by the 
MW. For equilibrium distribution of fractional occupan- 
cies fii, the anomalous term equals zero, but for occu- 
pancies fi2 > nie~ A2 > 1 ' Te induced by the MW resonance, 
Wa can lead to negative energy dissipation of the elec- 
tron system. In Eq. (|32[) . one can use the approximate 
expression for S' u , (q,u>i,i') given in Eq. (|2"3")l . Accord- 
ing to Eqs. (|23|) and (|32|) . the sign of —Wa coincides 
with the sign of ui2,i — (n' — n)uj c which can be nega- 
tive or positive depending on the magnetic field. Since 
^i,n;i',n' "C fc c < the contribution Wa is mostly 

exponentially small with the exception of magnetic fields 
where A;,;' - (n' - n) frw c < Ti tn -v,n'- 

The appearance of negative corrections to energy dis- 
sipation under the condition ni — n\e~ A2 < 1 / T ° > can be 
explained quite easily. The negative anomalous contri- 
bution (Wa > 0) corresponds to (n' — n)hjj c > A 2 .i. 
For narrow Landau levels, this means that scattering 
from the excited subband (/ = 2) to the ground sub- 
band (/' = 1) is accompanied by destruction of a ripplon, 
while the corresponding scattering back from the ground 
subband to the excited subband is accompanied by cre- 
ation of a ripplon. When ni = n,ie _A21 / Te , these two 
processes compensate each other in the expression for 
Wa- If n 2 > nie~ A2 - 1 / T ' , destruction of ripplons domi- 
nates, which leads to negative dissipation. In the oppo- 
site case, when (n — n) hio c < Aa ij creation of ripplons 
dominates, which results in additional positive dissipa- 
tion. It should be noted that the negative contribution 
to energy dissipation and negative momentum dissipation 



occur at the opposite sides of the point hu c = A2,i/m. 
Comparing Eq. (fTTj) with Eqs. (f2"5|) and (j2"5)l one can con- 
clude that the origin of this difference is the negative 
sign of the Doppler-shift correction in the ripplon excita- 
tion spectrum considered in the center-of-mass reference 



frame: E. 



(r) 



jq — ii^q ?iqV av . 

Consider now the energy loss rate of SEs due to elec- 
tron scattering at vapor atoms. In this case, the inter- 
action Hamiltonian is proportional to the density fluc- 
tuation operator of vapor atoms pk = 5Zk' a K'-K a K', 
where K = {K z , q} represents the momentum exchange 
between an electron and a scatterer. In terms of 
S$( ff ,w), the energy loss rate per an electron can be 
obtained as 



iK- 



l.V 



w = ~^y 

A 2 L\h Z - 



ni E ^K,K' 

1,1' K',K 
r( a ) • 



(33) 



where S^k,k' = £k< 



K 



Ej^, is the energy exchange 
at a collision. In order to obtain Wn and Wa, we shall 
firstly rewrite Eq. (|33l) trivially as a sum of two identical 
halves. Then, in the second half, the running indices 
K', K will be substituted as K'-K -> k', and K -> K, 
which changes the sign of the energy exchange, >^k.k' ~^ 
— . The next steps are the same as those resulting 
in Eq. (|29[) . Interchanging the running indices I <^ V in 

the second half, and using the basic property of sff, (q, u>) 
one can find 



W 



(y(»)) 



2A 2 L 2 h 



2ft E E 



^K,K' 



1,1' K',K 
xSi,i> (q,ui,i> ~ 
- - —A i IT 

ni — ni>e ' "t 



iK z z c 



N, 



(a) 
K' 



«tC,K') 



,(1/T-1/T e ) 



(34) 



This equation is more convenient for expansion in >?k,k' 
than Eq. (f33|) . 

Expanding Eq. (|34|) in >?k,k' , one can find again that 
W = Wn + Wa, where Wn and Wa have the same forms 
as that given in Eqs. (|30[) and (|3"2"|) . where x| r z ) should be 
substituted for 



~<a) m 
Xu ~ 



(V<- a ) 



2AL 2 TT e h 



(35) 

Using the condition K' 3> K, this equation can be sim- 
plified as 



(a) 



Xli 



(a) m e 
'° M 



En 



kUJr 



-Pl,l' 



(36) 



where 



" z 
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Expressions for C^ 1 and Bf^ convenient for numerical 
evaluations were given in Refs. and [ill 

The energy loss W transferred to vapor atoms and rip- 
plons is balanced by the energy taken from the MW field: 
W = (ni — n 2 ) A 2j ir, where r is the MW excitation rate 
defined by 



1 



n| 7 



(u — W2,l) +7 2 



(37) 



where 7 is the half-width of the MW resonance, and £Ir 
is the Rabi frequency proportional to the amplitude of 
the MW field. It is clear that negative contribution of 
Wa will be compensated by an increase in Wh due to 
electron heating. 

Some useful expressions for the SE energy relaxation 
rate obtained for arbitrary subband occupancies are 
given in the Appendix. It should be noted that nega- 
tive contributions to energy dissipation discussed above 
appear only for quasi-elastic iter-subband scattering. 
For SEs above superfluid 4 He, there are inelastic inter- 
subband scattering processes accompanying by simulta- 
neous emission of two short wave-length ripplon a 11 ' 17 . 
These processes cause strong additional energy relax- 
ation. Experiments of Refs. [If and were performed for 
SEs on the free surface of Fermi-liquid 3 He. For such a 
substrate, short wavelength capillary waves with q > 10 7 
are so heavily damped that even the existence of ripplons 
with such wave-numbers is doubtful. 



V. RESULTS AND DISCUSSIONS 

A. Vapor atom scattering regime 

Electron scattering at vapor atoms represents the most 
simple case for the magnetotransport theory, because the 
collision broadening of Landau levels of the same sub- 
band (Ti) is independent of the level number n. The 
same is obviously valid for the broadening of the gen- 
eralized factor Sim (9,w), which now can be denoted as 
T;.;. Additionally, the parameter defined in Eq. (j2"4")l has 



a very simple form A 



(a) 

Z,n;Z',n+m 



v^Pi,v (2n + 1 + m) 



which greatly simplifies evaluations. 

Consider electron temperature as a function of the 
magnetic field. It is defined by the energy balance equa- 
tion which contains the MW excitation rate r given in 
Eq. (|37p. In turn, r depends on the half- width of the 
MW resonance 7, which was studied theoretically with 
no magnetic field and under a parallel magnetic field 19 . 
If B is applied perpendicular to the surface, 7 should also 
have 1/B-oscillating terms, because inter-subband scat- 
tering increases when lS.i^jhLo c — > m. In our numerical 
evaluation, we shall use a qualitative extension of the 
result obtained for B = 0. According to this result, 7 
contains the contribution from intra-subband scattering 
722-11 and the contribution from inter-subband scatter- 
ing 72,1 = £"2->i/2. Under the magnetic field applied 
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FIG. 1. Electron temperature vs the magnetic field for three 
levels of MW irradiation estimated at B = 1 T: r = 10 5 s _1 
(1), 3 • 10 5 s _1 (2), and 5 • 10 5 s _1 (3). The solid curve was 
calculated for the model of 7(B) discussed in the text, the 
dashed curve represents the case r(B) = const. 



normally, electron scattering is enhanced by the factor 
hoj c / y/nTi 13 . Therefore, we can use an approximation 



722-11 



[Pa,2 - 2j?2,i] 



where r 2 ,i numerically is rather close to Y\ = 

K*J 2u c Vq^ /it. As for the oscillatory part entering 
72,1, we shall use the exact form of Eq. (|14|) . 

It should be noted that the oscillatory part of 7 is not 
large because p 2 ,i — 0.14. Still, it leads to some impor- 
tant consequences for electron temperature as a function 
of the magnetic field shown in Fig. [TJ Solid curves repre- 
sent results of numerical evaluations for the two-subband 
model taking into account oscillatory corrections to the 
MW resonance half-width 7, as described above. In this 
case, the electron temperature has small local minima at 
A2 i i/ftw c —¥ m due to oscillatory decrease in r a I/7. 
Three typical values of the Rabi frequency are chosen to 
provide MW excitation rate levels of 10 5 s _1 , 3 • 10 5 s _1 
and 5 • 10 5 s _1 at B = 1 T. For a model with a con- 
stant MW excitation rate r, which is applicable when in- 
homogeneous broadening dominates, the corresponding 
results are shown by dashed curves. At A 2 .i/fio; c w m 
these curves are nearly straight lines (without minima). 
For both models, the shape of curves describing the os- 
cillatory increase of electron temperature has asymmetry 
with regard to the point A2 i i/fia; c = m. This asymmetry 
is due to the negative correction of the anomalous term 
Wa leading to additional heating of the electron system 
at huj c > A2.1/771*. The asymmetry increases strongly 
with the MW excitation rate r and with m*(B). 

Electron heating increases with m* (lowering B), and, 
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FIG. 2. Electron temperature vs A2,i/ftaj c for the MW of 
higher resonant frequency, / = 130 GHz. Three levels of MW 
irradiation are the same as those described in the caption of 
Fig. [T] The all curves were calculated for the model of 7(B) 
discussed in the text. 



for the excitation rate r = 5 • 10 5 s" 1 at B = 1 T, the two- 
subband model fails at m* > 11. The applicability range 
of the two-subband model can be extended by using a 
stronger holding electric field which increases A2,].. In 
Fig. [2] electron temperature is shown as a function of 
the parameter A2,i/fiw c oc 1/B for a substantially higher 
MW frequency used in experiments on SEs 20 . For the 
solid curve, the two-subband model is applicable up to 
to* = 15. 

MW heating affects strongly the shape of conductivity 
oscillations because v c s depends on electron temperature. 
For example, the normal contribution to the effective col- 
lision frequency can be represented as 



(a) 

N, intra 



n lPl,l 



2^T e *f r M 



2T, 



coth U?r , (38) 



1A 



(a) 

N, inter 



vi a) h 2 -< 2 



Hij> can be approximated by a single term with to = to*, 
where m* depends on the magnetic field according to the 
above noted rule: to* = round (uji j>/u> c ). 

The anomalous contribution to the effective collision 
frequency has a different form 



» _ 



v\ = -- 



1-1/2 



i>u 



i%i - ripe 



-a^/tA Pl£_ 



r 2 

1 1:1' 



cot!) I ^ ) <!>,.,. 1^.1 + (-),.,- U; i 



(42) 



where functions (w c ) and 0;./' (co c ) are defined sim- 
ilar to Fim (oJc) and ffyj/ (uj c ) of Eqs. ((40]) and (|4T|) 
respectively, with the exception that their right sides 
contain the additional factor 2K(u)iji — muj c ) /Tiji orig- 
inated from Eq. (|23[) . Similar equations for SE energy 
relaxation rate are given in the Appendix. 



Comparing T e -dependencies of and u^' given in 



,(«) 



Eqs. (J38J), (J39J) and l|42|). we conclude that heating of the 
electron system reduces the normal contribution to the 
effective collision frequency. In contrast with this, the 
anomalous sign-changing correction , can be even en- 
hanced to a some extent with heating of SEs due to the 
factor coth (fkJ c /2T e ). Typical magnetoconductivity os- 
cillations of SEs calculated for the conditions of the ex- 
periment of Ref. [5j are shown in Fig. [3] Electron tem- 
perature calculated for these curves was shown in Fig. Q] 
by solid curves. At low MW excitation (r = 10 5 s _1 
at B = IT), magnetooscillations of a xx are just simple 
maxima centered at A2 j i/fiw c = to, which agrees with 
experimental observations. Between the commensura- 
bility conditions, a xx is suppressed, as compared to the 
dash-dot-dot line calculated for zero MW power. This 
suppression is due to 77.2 — > fi\ — ¥ 1/2 and weaker scat- 
tering at the excited subband. The increase in the decay 
rate ^2^1 at A2,i/fto; c — > to leads to a sharp decrease in 
H2, which restores a xx values obtained without the MW 
field. This restoration is not complete if T e > T, because 



,(«) 

~N,intr 



a decreases with heating, as discussed above. The 



'A tl UJ •s—^ I A l T \Pll> 

" r £ + ni ' e ) r~ x At hi S her MW excitation 

2y/TTl e ., \ ' i-lj' at H - 1 Tl flip oliano nf n 



l>l 



coth ( — ) Fiji (w c ) + Hi tV (w c ) 



where the new functions 



Fi,i> (cj c ) = 22 ex P 



Him (w c ) = ^2 TOexp 



ft 2 {uiji - muj c y 
f 2 



H 2 (uij> - muJcY 
Y 2 

1 Id' 



(39) 



(40) 



(41) 



v n \ inter ^ s ver y small under these conditions. 

= S-IOS" 1 and 5- 10 s s- 1 
at B = 1 T), the shape of conductivity oscillations is af- 
fected strongly by the anomalous term ir^ leading to lo- 
cal minima at A2,i/fiw c > to. The increases with r 
because of two reasons. The first reason is the increase in 



An = ri2— n\e 



-A 2 ,i/T e 



at higher MW excitation shown in 



defined for I > I' are independent of T e . For narrow 
Landau levels (Fiji <C hu) c ), the series defining Fiji or 



Fig. 0J The second reason is electron heating due to decay 
of the excited SE state which increases coth (huj c /2T e ) of 
Eq. (|4"2")l . Further shape evolution is shown in Fig. [5] 
for lager to* and / = 130 GHz, where, according to 
Fig. [3J electron temperature can take a value of about 
2 K. As expected, under these conditions the anomalous 
contribution strongly increases. Near the commensura- 
bility conditions l\i.\jtvjj c — > to, the shape of conductiv- 
ity oscillations becomes similar to that observed for the 
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FIG. 3. Magnetoconductivity vs A 2 ,i/ftu c oc 1/B for the MW 
of / = 79 GHz. The dashed-dot-dot line was calculated with 
no MW irradiation. Dotted, dashed and solid curves represent 
three levels of MW irradiation (r) given in the caption of 

Fig.m 



FIG. 5. Magnetoconductivity vs A.2,i/hw c for the MW field of 
higher resonant frequency, / = 130 GHz. The dashed-dot-dot 
line was calculated with no MW irradiation. Dotted, dashed 
and solid curves represent three levels of MW irradiation (r) 
the same as those in Fig. [2] 




0.59 



FIG. 4. Deviation of subband occupancies from the equilib- 
rium distribution, An = n% — fiie _A2 ' 1 ' /Tc , vs the parameter 
A2,i/&j c oc 1/B. Dotted, dashed and solid curves represent 
three levels of MW irradiation (r) given in the caption of 

Fig.m 

electron-ripplon scattering regime 5 at T = 0.2 K. It is 
important that the results given in Fig. [5] are still ob- 
tained in the validity range of the two-subband model. 

It is instructive to compare the peak broadening of dif- 
ferent quantities shown in Fig. [5] The broadening of the 
decay rate i>2^i coincides with T2 l i, which is an aver- 
age of T 2 and Ti. In contrast, such quantities as |z/a|, 
fii — 1/2, and T e — T rise in a much broader magnetic 
field range having nearly the same width which does 



FIG. 6. Line shapes of i>2^i (dashed), va (solid), T e — T 
(short-dotted), and fii — 1/2 (short-dashed) as functions of 
B near the commensurability point with m = 5, under the 
conditions: T = 0.4 K, and r = 3 • 10 5 s -1 . 

not represent the broadening of Landau levels directly. 
We shall use this similarity in the line widths of T e — T 
and fix — 1/2 later, considering electron heating for the 
electron-ripplon scattering regime. 

B. Electron-ripplon scattering regime 

For electron-ripplon scattering, the anomalous (sign- 
changing) contribution to the effective collision frequency 
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FIG. 7. Evolution of the a xx (B) line shape near commensu- 
rability points with the gradual increase in m* at T = 0.2 K 
and / = 79 GHz. 



is induced by the MW resonance absolutely in the same 
way, as it is for electron scattering at vapor atoms. In the 
case of liquid 3 He, electron-ripplon scattering dominates 
at low temperatures T < 0.2, where the half-width of the 
MW resonance 7 is substantially reduced. According to 
Eq. (137)) , at the same amplitude of the MW field, this 
decrease in 7 leads to a strong increase in the MW ex- 
citation rate r at the resonance u — cj 2 ,i, which greatly 
magnifies v^. 

Unfortunately, the electron-ripplon scattering regime 
is much more difficult for the analysis of the effect of 
electron heating than the vapor atom scattering regime 
because of different reasons. First, the electron-ripplon 
coupling Uq has a very complicated form 1 *: 



— - Ki (qz) 
qz 



eE± 



dV e 



(0) 



3z e 



where K\ (x) is the modified Bcssel function of the sec- 
ond kind, and Ve°^ (z) is the electron potential energy 
over a flat surface. Therefore, it is impossible to obtain 
simple analytical equations for the energy loss function 
W(T e ). Moreover, if 3 He is used as the liquid substrate, 
there might be contributions from other mechanisms of 
energy relaxation, which by now have no strict theoreti- 
cal descriptions. 

As indicated above, heating of SEs only increases the 
importance of the anomalous contribution to the effective 
collision frequency. For electron-ripplon scattering, it fol- 
lows directly from Eqs. (|25 [) -(|27 p . Therefore, in order to 
prove the possibility of existence of zero resistance states 
due to non-equilibrium filling of the excited subband, it 
is sufficient to show that negative a xx can appear even 
without electron heating. At T e = T = 0.2 K, the MW 
field amplitude, which gave r = 10 5 s _1 (at B = IT) 



for the vapor atom scattering regime shown in Fig. [3] 
(dotted curve), now gives r = 2 • 10 6 s _1 , because the 
MW resonance line width 27 ~ 0.3 GHz due to inhomo- 
geneous broadening^. This excitation rate is very high, 
because it leads to a xx < already at m* = 4. For 
presentation of Fig. [7J we had chosen a two-times lower 
excitation rate r = 10 6 s~ 1 independent of the magnetic 
field. This figure shows the evolution of the line shape of 
conductivity oscillations with the gradual increase in the 
integer parameter m*. It is quite convincing that even 
without heating of SEs the anomalous contribution to 
the effective collision frequency increases strongly with 
m* , and the conductivity curve corresponding to m* = 6 
enters the negative conductivity regime in the vicinity of 
the minimum. This is in accordance with experimental 
observations reported for the high magnetic field range 
(m* < 10). 

Maxima and minima of ua(B) have the same ampli- 
tude. Without heating of SEs, amplitudes of conductiv- 
ity maxima obtained here are larger than amplitudes of 
minima, because the normal contribution is^ increases 
at A.2 t i/ hu) c — !• m due to oscillations of subband occu- 
pancies. Experimental curves^ show that at strong MW 
power and large to* amplitudes of minima are larger. 
This could be an indication of electron heating, because 
t^j ~ 1/T e . To analyze the effect of heating of SEs on 
conductivity oscillations, we shall model electron temper- 
ature oscillations using similarities in the line shapes of 
T e — T and Ui — 1/2 shown in Fig. [6] In particular, we 
assume that an electron temperature peak is described 
by T e (B) = T + 2 (AT e ) max [nx(B, T) - 1/2], where the 
maximum elevation (AT e ) max depends of to*. We disre- 
gard the asymmetry of the peak induced by W\ because 
it does not lead to a substantial change in final results. 
The results of such a model treatment of the heating 
effect are shown in Fig. [8] They indicate that even mod- 
erate heating of SEs affects strongly the shape of magne- 
tooscillations, making amplitudes of minima larger than 
amplitudes of maxima (dotted curve) in accordance with 
experimental data. 

In Fig. El we had chosen the excitation rate r = 
10 6 s _1 , so that the initial curve (solid) calculated for 
T e = T have a small minima with a xx > 0. Then, we 
found that heating with (AT) max = 0.1K strongly re- 
duces conductivity extremes due to ^ ~ 1/T e , and mod- 
erate heating with (AT) max = 0.5 K leads to a minimum 
with a xx < 0. Therefore, decay heating of electrons, 
which occurs in the vicinity of the commensurability con- 
ditions, helps to obtain zero resistance states. For exam- 
ple, within the validity range of the two-subband model, 
electron temperature peaks of about 2 K can reduce j/n 
by an order of magnitude. Still, heating alone cannot 
make a xx < 0. It is only the anomalous contribution i/a 
which eventually leads to negative conductivity and zero- 
resistance states. Without i/a, a conductivity dip would 
be an even function of the parameter u>2,i — rn*to c with 
&xx > 0. The existence of a magnetoconductivity max- 
ima at the opposite side of the point W2,i ~ m*uj c = in 
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bic correction increases Yi >n .i' >n i by about 1.3, if we as- 
sume x q ~ 1. If we take into account that the integrand 
of Eq. (2U has a maximum at x q ~ m* + 2, the broad- 
ening increases approximately two times. Therefore, a 
qualitative analysis indicates that the many-electron ef- 
fect becomes more important in the low magnetic field 
range where it increases the width of conductivity oscil- 
lations and reduces amplitudes of maxima and minima, 
which also agrees with experimental observations. Decay 
heating increases the Coulombic correction to the broad- 
ening of magnetooscillations. Still, a strict description 
of Coulombic effects on magnetoconductivity oscillations 
requires a more careful study. 



VI. CONCLUSION 



FIG. 8. Evolution of the a xx (B) line shape near the com- 
mensurability point m* = 4 with the increase in (AT e ) max 
at T — 0.2 K and / = 79 GHz: (AT e ) max = (solid), 0.1K 
(dashed), and 0.5 K (dotted). 



experimental curves, which demonstrate vanishing mag- 
netoconductivity^, is an additional evidence for a sign- 
changing correction convincing that ZRS are realized at 
the vanishing points. 

It should be noted that at T — 0.2 K, one-ripplon scat- 
tering processes are not sufficient to prevent strong heat- 
ing of the electron system at the commensurability con- 
ditions. In particular, for r = 5- 10 5 s" 1 , estimation gives 
(T e ) max ~ 3K at m* = 4. The model treatment of the 
heating effect discussed here allows to draw conclusions 
about actual role of the electron heating in experiments 
with SEs 5 . For example, the firm conductivity maxi- 
mum (without a minima) observed for radiation power 
P of -25 dB at m* = 4 surely indicates that electron 
heating is small or moderate under these conditions, and 
there is an additional mechanism of energy relaxation 
at low ambient temperatures. We speculate, that the 
magnetopolaronic effect and electron coupling with bulk 
quasi-particles, giving a very small correction to the mo- 
mentum relaxation rate under experimental conditions, 
can contribute substantially to the energy relaxation rate 
reducing electron temperature. 

Experiments^ 5 - are conducted for low surface electron 
densities n s of about 10 6 cm -2 . Nevertheless, electron- 
electron interaction affects noticeably experimental data. 
According to Ref. Hi], under magnetic field an elec- 
tron moves in a quasi-uniform electric field of other 
electrons Ef of fluctuational origin. Its average value 
E^p ~ i^/Tln^J 4 increases strongly with electron tem- 
perature and density. The fluctuational electric field 
increases the broadening of the DS F 11 ' 16 Ti t n;i',n' - > 

\l T ln;i>,n> + where T c - y/2eEfl B oc 1/y/B. 

Thus, at T e = 0.2 K, and n s = 0.9-10 6 cm" 2 , the Coulom- 



In summary, we have developed the theory of magne- 
toconductivity oscillations in a multi-subband 2D elec- 
tron system under MW irradiation of a resonant fre- 
quency. We have shown that besides the quite obvious 
1 /B-modulation of conductivity, the non-equilibrium fill- 
ing of the excited subband induced by the MW resonance 
leads also to sign-changing corrections to the effective col- 
lision frequency due to usual inter-subband scattering. 
As the MW power goes up, the corresponding increase 
in the amplitude of these sign-changing corrections can 
result in the negative linear response conductivity and 
zero-resistance states. 

Our theory is based on the self-consistent Born approx- 
imations, and it is presented in a general way applicable 
for any quasi-elastic scattering mechanism. As particu- 
lar examples, we have considered two kinds of scatter- 
ers which are typical for the electron system formed on 
the free surface of liquid helium: helium vapor atoms 
and capillary wave quanta (ripplons) . In the vapor atom 
scattering regime, we found a strong 1/B-modulation 
of the electron temperature, which increases sharply in 
the vicinity of commensurability conditions. This de- 
cay heating is shown to enhance the effect of the sign- 
changing terms in the longitudinal conductivity <j xx . The 
evolution of the line-shape of conductivity oscillations 
with an increase of the MW field amplitude is studied, 
taking into account heating of surface electrons. 

For the electron-ripplon scattering regime, we have 
shown that magnetooscillations of large amplitude and 
the negative linear response conductivity of SEs can eas- 
ily appear under moderate MW excitation even for cold 
SEs. The evolution of the line-shape of a xx extremes 
caused by an increase in the electron temperature is stud- 
ied using a model treatment. We believe that theoretical 
results presented in this work explain all major features 
of MW-resonance-induced magnetooscillations observed 
in the system of SEs on liquid helium, and support the 
suggestion^ that novel zero-resistance states are realized 
in such a system. 
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Appendix A: Energy relaxation rate 

Here we give final expressions for the energy relax- 
ation rate of SEs due to scattering with vapor atoms. 

The normal 9^ and anomalous energy relaxation 
rates are defined by the following relationships: Wist = 



^SLro + Winter ' where 



(T e -T)4 a) , and W A 



In turn, v 



N 



-(a) 
N Antra 



n l/2 MT 1^ ~ 



1:1 



u M + _p M coth^— 



(Al) 



-(a) _ VQ a) m e frw c E R ni + n v e A '.''/ Te 
"Winter- n l/2 MT 1^ T 

e 1>V M 



v 



uij'Fiji (LJc) + -=-pi,v 
Er 



Fiji (w c ) coth 



fvjj c 
2T, 



+ Him (w c ) 



, (A2) 



I 

functions Fn> and H\ y were given in Eqs. (|40[) and (1411) . The anomalous energy relaxation rate can be repre- 

sented as 

I 



_(„) _ 2v { a 'm e huj c E R ^ ("/ -ni>e A '.i'/ T =) 



ttV2M 



r 2 



UlV®n> (UJ C ) + r=-ptl> 
Er 



<!>/./' U',-)cotl) ( ) + (-'•' 



(A3) 



I 

where &ij> andG/,;' are the same as those of Eq. (|4^|) . For tions were used for obtaining electron temperature as a 
equilibrium subband occupancies, i>a = 0. These equa- function of the magnetic field under the MW resonance. 
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